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'. Abstract 

O 

*q | We obtain an approximate global stationary and axisymmetric so- 

O^ 1 lution of Einstein's equations which can be considered as a simple star 

5h . model: a self-gravitating perfect fluid ball with constant mass den- 

sity rotating in rigid motion. Using the post-Minkowskian formalism 
(weak-field approximation) and considering rotation as a perturbation 
^ . (slow-rotation approximation), we find approximate interior and ex- 

terior (asymptotically flat) solutions to this problem in harmonic and 
quo-harmonic coordinates. In both cases, interior and exterior solu- 
tions are matched, in the sense of Lichnerowicz, on the surface of zero 
pressure to obtain a global solution. The resulting metric depends on 
three arbitrary constants: mass density, rotational velocity and the 
star radius at the non-rotation limit. The mass, angular momentum, 
quadrupole moment and other constants of the exterior metric are de- 
termined by these three parameters. It is easy to show that this type 
of fluid cannot be a source of the Kerr metric. 

keywords: relativistic astrophysics, two post-minkowskian approxima- 
tion, harmonic coordinates, rotating stars 



1 



1 Introduction 



Rotating stars have been an active field of research in General Relativity in 
the last sixty years (see [H El El IU E] and references there in). Nevertheless, 
until the discovery of neutron stars, this research was mainly theoretical since 
the relativistic effects of rotation are negligible for normal stars. The high 
angular velocit and small radius of neutron stars have stimulated new efforts 
to understand the gravitational fields outside and inside a rotating star in a 
non-Newtonian framework. 

Only a few exact solutions of the Einstein equations can describe the 
gravitational field inside a rotating star, i. e., the gravitational field of an 
axisymmetric rotating perfect fluid (see [0] for a brief review and [HE])- So 
far it has not been possible to match any of them to an asymptotically flat 
vacuum solution, the gravitational field outside the star, to obtain a global 
star model similar to those provided by Newtonian theory. This lack of ex- 
amples does not mean the absence of theoretical results to this problem even 
in full nonlinear General Relativity. For instance, axisymmetric matching 
conditions have been largely revised and used to demonstrate the uniqueness 
of the exterior gravitational field to an isolated body |9J. Nevertheless the 
most relevant outcomes result from approximation methods and numerical 
computation. Both techniques have been applied to obtain useful descrip- 
tions of the gravitational field of a star obeying different equations of state 
(relativistic polytropes, hadronic fluids, etc. . . ). They have also been used to 
study other related problems such as stability, normal modes of vibration and 
generation of gravitational waves (all these topics are extensively reviewed 
in0). 

We are interested in analytic approximation methods. Obviously, mod- 
els that can be handled with these methods have to be simpler than those 
amenable to numerical methods. Analytic approaches provide a better un- 
derstanding of the basic features of relativistic rotating stars, so they be of 
some help to study more interesting and realistic models which deserve more 
powerful methods. 

On these bases we introduce a new analytic approximation scheme and 
apply it to search for an approximate global solution to the gravitational 
field of a fluid with a very simple equation of state. We aim to establish an 
approximate solution of the Einstein equations which describes the gravita- 
tional field inside a ball of perfect fluid with constant mass density rotating 
in a rigid motion, and to match it, on the zero pressure surface, to an asymp- 
totically flat approximate solution of the vacuum Einstein equations. We 
solve this problem in harmonic coordinates and quo-harmonic coordinates 
(harmonic coordinates in the sense of the quotient metric). In both cases, 
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we require the metric to be of class C 1 on the matching surface, that is the 
metric should satisfy the Lichnerowicz matching conditions |10j . 

Two classes of analytic approximation methods are described in the lit- 
erature. First, the post-Newtonian approximation scheme, which computes 
perturbations to well-known Newtonian solutions such as the McLaurin ellip- 
soids [TTJ [T2] ; second, the slow rotation approach, which calculates perturba- 
tions to spherically symmetric solutions of the Einstein equations |13j . The 
approximation scheme we propose is an intermediate way: we implement a 
slow rotation approximation on a post-Minkowskian algorithm (a first at- 
tempt can be found in [2]). We introduce two dimensionless parameters. 
One, A, measures the strength of the gravitational field, the other, Q, mea- 
sures the speed of rotation of the fluid. We also make some assumptions 
about how all the constants that appear in the post-Minkowskian metric de- 
pend on these two parameters. So we obtain an expansion of the metric in a 
double power series of A and Q. Obviously, if there is no rotation ($1 = 0), we 
are faced to the post-Minkowskian perturbation to the Newtonian gravita- 
tional field of a spherically symmetric mass distribution. On the other hand, 
Newtonian deformation of the source due to rotation is included in first order 
A terms up to some order in the rotation parameter. In other words, there 
is no exact Newtonian rotational model underlying our approach but only 
an approximate one. Nevertheless, we use these linear terms to compute 
non-linear contributions of rotation to the gravitational field. 

We only evaluate metric terms of order less than or equal to Q 3 and no 
greater than A 2 . However, since the algorithm is implemented by an algebraic 
computational program, our results can easily be enhanced, if so desired, by 
going farther in the approximation scheme. The values taken by the two 
control parameters for a definite system tell us which terms in the double 
expansion of the gravitational field are relevant to that particular problem. 

In section 2 we list the properties we impose on the space-time and the 
energy-momentum tensor and we define the notation. Moreover we briefly 
review the post-Minkowskian approximation we use in the following sections. 

In section 3 we obtain the first order post-Minkowskian metric up to 
the second order in the rotation parameter. We solve the linearized Ein- 
stein equations outside and inside the mass distribution. We discuss the 
dependence of the free constants in the linear solution on A and Q. We use 
that information to cancel all terms in the solution of order higher than O 3 . 
Finally we match those simplified exterior and interior solutions on the ap- 
propriate surface up to this order of approximation. Consequently we derive 
expressions for all the surviving free constants of the exterior and interior 
gravitational fields as functions of the parameters of the source, its mass 
density, mean radius and angular velocity. A global system of harmonic 
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coordinates is assumed. 

In section 4 we work up the second-order approximation. 

The role of frames of reference in General Relativity [15] leads us to focus 
our attention on the quo-harmonic coordinates. In section 5 we show that 
our scheme can also be successfully applied by imposing on the metric this 
unusual but interesting coordinate condition. 

2 Preliminaries: notation and other topics 
2.1 Metric 

We require space-time to be a stationary and axisymmetric Riemannian man- 
ifold admitting a global system of spherical-like coordinates {£, r, 9, y?} which 
verifies the following properties: 

1. Coordinates are adapted to the space-time symmetry, £ = d t and £ = 
d v are respectively the timelike and spacelike Killing vectors, so that 
the metric components do not depend on coordinates t or if. 

2. Coordinates {r, 8} parametrize two dimensional surfaces orthogonal to 
the orbits of the symmetry group, that is the metric tensor has Papa- 
petrou structure, 

g = 7« uj 1 ® u* + 7 t ^(V ® ^ + uj^p® u l ) + 7 W u v ® 
+ 7 rr u r ® uj r + 7 re (a; r <g> uj + uj% u r ) + j 0e u e ® u e , (1) 

where oj* = dt, uf = dr, uj e = r dd, uj v = r sin 9 dip is the Euclidean 
orthonormal cobasis associated to these coordinates. 

3. Coordinates {t, x = r sin 6* cosy?, y = rsin^siny?, z = cos^} associ- 
ated with the spherical-like coordinates are Cartesian coordinates at 
spacelike infinity, that is the metric in these coordinates tends to the 
Minkowsky metric in standard Cartesian coordinates for large values 
of the coordinate r. 

All these properties are compatible with the two classes of coordinates 
we use in this paper, mainly harmonic and quo-harmonic coordinates. Time 
coordinate t is always harmonic under assumptions 1 and 2, then we have 
only to check that spatial coordinates {x, y, z) verify the coordinate condition 
we choose. 
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2.2 Energy- momentum tensor 



We assume that the source of the gravitational field is a perfect fluid, 

T={^ + p)u®u + pg, (2) 

whose density /x is a constant and whose pressure p depends only on r and 9 
coordinates. We also assume that the fluid has no convective motion, so its 
velocity u lies on the plane spanned by the two Killing vectors, 

u = V(£ + a;C), (3) 

where 

ip = — 0y u + 2uo j t(p r sin 6 + uo 2 7 W r 2 sin 2 6j 2 (4) 

is a normalization factor. 

Let us consider Euler equations for the fluid (or, what is equivalent, the 
energy-momentum tensor conservation law) to gain information about the 
matching surface. Assuming the fluid rotates rigidly as we do, u = constant, 
the equations read 



d a p = (fi + p)d a In ip (a,b, . . . = r ,6) . (5) 

Since \x depends neither on r nor on 9, we can integrate them to give a very 
simple expression for the pressure in terms of the normalization factor 

which in turn leads to the following implicit equation for the matching surface 

p = ^ V = ^e, (7) 

here ip-z is an arbitrary constant. 

These last two equations, © and ((7|), play an important role in our 
scheme. We use them to derive approximate expressions for the pressure and 
the matching surface in a coherent way with the expansion of the metric we 
propose. 

2.3 Einstein's equations: Post-Minkowskian approxi- 
mation 

The Einstein equations for any energy-momentum tensor read 

Rap = 87T (T al3 - ^g a pT^J = 87Tt a/3 . (8) 
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The deviation of the metric tensor from the Minkowski metric, 



h a p = 9a(3 — Va 



(9) 



in a system of Cartesian-like coordinates, where (7} a p) = diag(— 1, 1, 1, 1), 
provides the natural object to work out a post-Minkowskian approximation. 
If we substitute the above expression for the metric into equation (jHJ) and 
neglect all the quadratic and higher terms in h a p, we obtain the well known 
linear approximation to the Einstein equations p2] , but if we do not neglect 
these terms we can use equation (jSJ) to generate a set of successive approx- 
imations to the exact solution g a p. This approximation method is usually 
implemented in harmonic coordinates, 



(10) 



although another type of coordinates could be used (for instance, those we 
shall use in section 5, i. e., quo-harmonic coordinates). 

Using the metric deviation, the Ricci tensor can be written as 



R 



a/3 



1 . l n 

--nh af3 + -d Q 



d p h 



ha 



1 , 



+ 5% 



d p h 



b pa 



a/3 ) 



ir 



where □ = r} a/3 d a dp , h = r] al3 h a/ 3 , d p = f] p,J d a and N a p collects quadratic 
and higher terms in h a p. Analogously the harmonic condition (TlOT) takes the 
form 

(12) 
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v P a9 Xfi r p = d p h pa - ~r) pa h)+H a = o, 



Xp - \ --fm 2' 

where H a also takes account of non linear terms. These equations reduce to 



Ah a/3 = 2 {N af3 - 8irt al3 - d( a Hfl) = S, 



a/3 , 



d 3 ( h ja - -i] ja h ) = -H a 



(13) 



when we look for stationary metrics, doh a p = 0. 

Let us remark that something must be said about how the energy-momen- 
tum tensor depends on the metric in order to obtain a true perturbation 
scheme from the above equations. Such information is given by formulas (ED, 
(j3J), (BJ and (jHJ) for the type of fluids we are considering in this paper if we 
identify the Cartesian-like coordinates {t, x, y, z} of our metric with the 
coordinates {x a } used in the post-Minkowskian algorithm. 
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3 First order metric in harmonic coordinates 



3.1 Linear exterior solution 

In the absence of matter, t a p = 0, equations ( fl3l) becomes linear if all the 
quadratic and higher order terms in the metric deviation h a p are neglected, 
that is N a/ 3 = H a = 0, 

Ah a/ 3 = , 

d k (h kfl -^hri^) = 0. (14) 

The general axisymmetric solution of these homogeneous equations which is 
compatible with the Papapetrou structure we have assumed for the metric 
[see formula (0Q)] and has a regular behavior at infinity is well known p^8l fT9] . 
It can be written as follows: 

00 M oo j 

fc = 2£^(r,, + D,,) + 2i;^z„ 

n=0 ' n=l ' 

+ T.^^ + ^^- 1 F n , (15) 

n=l ' n=2 ' 

where 

T n = P n (cos 6) iJ ® (J (n > 0) , 
D n = P n (cos 9) 5ijdx l ® dx j (n > 0) , 

Z n =P n 1 (cos^)(a;' ®u* + u ip ®u> t ) (n > 1) , (16) 
are spherical harmonic tensors and 

E n = ~nD n + ^nHl + H 1 n (n > 1) , 

F n = \n{n-l)Hl + {n-l)Hl-\Hl (n > 2) (17) 

are two suitable combinations of and these other three spherical harmonic 
tensors 

Hi = P n (cos 9) (5ij - Se^dx*® dx j (n > 0) , 
i/^ = P^(cos 9) (kitj + kje^dx 1 ® dx j (n > 1) , 
//; = P n 2 (cos 0) {kikj - rriim^dx*® dx j (n > 2) (18) 

ki, €i and mi stand for Euclidean unit vectors of standard cylindrical coor- 
dinates, dp = kidx 1 , dz = eidx 1 , pdip = m^dx 1 (this is the set of spherical 
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harmonic tensors we use to write covariant tensors of rank-2 in this paper); 
M n and J n are the multipole moments of Thorne [18J or Geroch-Hansen [20] 
(both definitions are equivalent [21]) and A n and B n are other constants that 
can be related to "stress moments" of the source [22J . These last two sets of 
constants unlike M n and J n are not intrinsic, but they are necessary to solve 
the Lichnerowicz matching problem. 

F 2 has spherical symmetry (therefore it must be included in the spherical 
symmetric linear solution besides the mass monopole term). This is easy to 
verify from the following alternative expression (indices are raised with the 
Euclidean metric) 

if = H° 2 ij + Hi ij - X -H\ ij = 6 lj - 3nV" , (19) 

where n l is the Euclidean unit radial vector of standard spherical coordinates. 

We want the metric to have equatorial symmetry, so we set M 2n +i = , 
J 2n = , A 2n+ i = and B 2n+ i = in equation ([TBI) . 

3.2 Multipole moments and control parameters 

This section is the key to understandig the current analysis. We expect all 
those constants of the exterior metric (1151) to be functions of the parameters 
included in the energy-momentum tensor of the fluid; i. e., the mass density 
the angular velocity to and the value ips of the normalization factor on the 
zero pressure surface. Before making any assumption on the dependence of 
those exterior constants on these parameters, let us introduce another set of 
related parameters which will be more suitable to this end. 

First let us denote by r the typical size or mean radius of the source, it 
can also be viewed as the radius of the sphere towards which the matching 
surface would tend at the non-rotating limit. Therefore m = ^[ir^ may 
provide a good estimate of the total mass of the source and the dimensionless 
parameter A = m/ro (the ratio of Schwarzschild's radius to the mean radius 
of the source) may give a measure of the strength of the gravitational field. 
We shall assume that this parameter is small and use it to control the post- 
Minkowskian expansion of the metric, in the sense that none of the terms of 
the linear exterior solution can be of order A 2 or greater. 

Now let us take a look to the implicit equation of the matching surface. 
We expect that expression ([7]) defines a slightly deformed sphere, r ps r , 
at the lowest order of approximation. Nevertheless, this is not compatible 
with the post-Minkowskian approach, because it imposes at this level 7# ~ 
— 1 + 0(A), "/tip ~ 0(A) and 7 W ~ 1 + 0(A) what leads to a cylinder, 
p = constant, not to a sphere, unless we assume that u 2 is at least of order 
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A. That suggests that we should introduce another dimensionless parameter 
Q, such that Q 2 = TqU 2 /X, instead of the apparently more physical parameter 
tq to. This new parameter can be seen as an estimate of the ratio of centrifugal 
to gravitational forces or rotational to gravitational energies on the surface 
of the source [23]. At all extends Q is a parameter that takes account of the 
deviation of the solution from spherical symmetry due to rotation. 

Once we have decided the parameters we are going to use in our approach, 
we have to elucidate how the constants in (1151) depend on them. Our proposal 
can be summarized in the following set of substitutions (we use the same 
notation for old and new constants): 

M n -> mr™M n = \r™ +1 M n -> Afi n r" +1 M n , 
J n - mur n +l J n = A 3 / 2 fir " +1 J n - A 3 / 2 O n r£ +1 J n , 

B n -> \r% +1 B n -> \n n - 2 rZ +1 B n . (20) 

The first step drives multipole moments to a dimensionless form. An extra 
factor Q is included in all J n multipole moments due to their clear dynamical 
character, they should vanish if there is no rotational motion of the fluid. 
This feature has also been taken into account in the second step, since all 
multipole moments except M and B 2 (as it was pointed out at the end of the 
preceding section) should also vanish if the solution has spherical symmetry. 
The leading power of Q we assign to each constant in this second set of 
substitutions reflects just the order of the spherical harmonic polynomial to 
which it is related. This is partly explained by the structure of the mass 
multipole moments of the MacLaurin ellipsoids: the leading term in the 
expansion in powers of Q of these quantities shows the same dependence on 
this parameter [14J. 

Substituting ( 1201) into expression ( fT5i) and neglecting terms of order equal 
or greater than f2 4 (slow rotation) we arrive at the following approximate 
expression for the exterior metric 

g ext ^ -1+2A-S T + 1+2A-^ D + A-|F 2 
V V J V V J T 

+ 2\n 2 -^ (T 2 + D 2 ) + \n 2 ^ E 2 + Afi 2 ^ F 4 

+ 2A 3 / 2 ^4 Zx + 2A 3 / 2 fi 3 4 ^3 , (21) 

where r] = r/r is a dimensionless variable. Let us remark that up to this 
order of approximation M 2 , J3, A 2 and B4 are pure numbers but Mq, J\ and 
B 2 can be linear functions of Q 2 . 
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3.3 Matching surface and energy- momentum tensor 

Let us assume the interior metric has the same dependence on A than the 
exterior metric (EH) has, 



7« « -1 + Xfu , lt v ~ A 3/2 /^ , 7 W »1 + A/ w , 

7 rr « 1 + Xfrr , 7r0 ~ Xf r9 , 700 ~ 1 + ( 22 ) 

(we shall not use labels to distinguish between exterior or interior metrics 
not even between first or second order approximate metrics whenever it can 
be clearly understood to which of them we are referring). Thus the normal- 
ization factor (jlj) up to first order in A reads, 

^^l + ^\(f tt + n 2 V 2 sm 2 9) . (23) 

As far we are assuming that the metric components are continuous on the 
matching surface, we use their exterior expressions given by (12T]) to make up 
(17j) into a true equation for such a surface. Then we search for a parametric 
form of the matching surface up to zeroth order in A by making the following 
assumption, 

r « r (l + att 2 P 2 (cos9)) . (24) 
A simple calculation leads to 

"-Sbfa-D- (25) 

We also obtain a similar expression for in terms of the exterior constants, 

w 1 + A (m + ^ . (26) 

To complete the linear interior equations we need an approximate expres- 
sion of the energy-momentum tensor of the fluid. First of all, let us notice 
that the density fi = 3A/(47rrQ) is a quantity of order A. Thus taking into 
account ( 1231) and ( 1261 . it is easy to check that the pressure (jHJ) is of order A 2 . 
Therefore the energy-momentum tensor (j2J) contributes to the right side of 
the Einstein equations by means of 

A A 3 / 2 f2 
Stt t « 3- (T + D ) + Q^^v Zi , (27) 

if terms of order equal or higher than A 2 are disregarded. 
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3.4 Linear interior solution 



Let us consider the following system of linear differential equations 

A/i a/3 = -16irt a f3 , 

d k (h kfi -^hrj kfl ) = 0, (28) 

where t is given by (12 7p . The general solution of these equations which is 
regular at the origin of the coordinate system r = and has the Papapetrou 
structure (OQ) is 

oo oo 

h = h inh + £ m„ r n (T n + D n ) + J2 3n r n Z n 

n=0 n=l 

oo oo 

+ a D + b H + ]T a n r n E* n + £ b n r n F\ , (29) 

n=l n=2 

where 

h inh = -Xr, 2 (T + D ) - ^A 3/2 fir/ 3 Z x (30) 

is a solution of the inhomogeneous system. We have also introduced two new 
sets of spherical harmonic tensors, 

E\ = ^{n + 1) D n - + l)H n + H l n {n > 1) , 

F* B = i(n + l)(n + 2)H n -(n + 2)Hi-iHj (n > 2) , (31) 

which seem to be suited to the interior problem better than those we used 
to write the linear exterior metric. 

The integration constants that appear in expression (I2"9l should have 
some dependence on the physical parameters A, Q and r similar to that 
we propose for the constants of the exterior metric (there we used capital 
letters to label the constants, here we use small case letters). We assign to 
m n , j n and a n the same dependence on A and Q that we assumed for their 
corresponding capitals. This rule must be slightly modified for the constants 
b n , in the sense that they differ from their capitals in a factor Q 2 (F* 2 is not 
a spherical symmetric tensor). There are two exceptions, ao and bo, because 
they have not counterpart in the exterior metric, we assume that both are 
linear in A and zeroth order in Q (in fact b must be of second order in Q 
since H has not spherical symmetry). Finally, we must include a factor r^™ 
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in all constants with an index n as they are associated to r n terms of the 
interior metric. It can be summarized as follows: 

m -> Am , m 2 -> Afi -=- , ji -> A 7 fi— , j 3 -> A 7 Q , 

?o r rg 

a — > Aa , 6 ^A6 , a 2 — > Afi 2 ^| , 6 2 ^Afi 2 -|. (32) 

Constants we have not mentioned above whether depend on higher powers of 
Q or are incompatible with the equatorial symmetry requirement, so they can 
be ignored. As before, these new dimensionless constants are pure numbers, 
m 2, hi CL2 and b 2 , or linear functions of Q 2 , m , ji, a and b . 

Finally, substituting these values for the constants into equation (129|) . we 
arrive to the following approximate expression for the interior metric: 

g int « (-1 + Am - Xrf) T + (1 + Am + Xa - X V 2 ) D + A6 H Q 
+ \n 2 m 2 r] 2 (T 2 + D 2 ) + Xn 2 a 2V 2 E* 2 + Xn 2 b 2 r] 2 F* 2 

+ X 3/2 £lr] (j x - Z, + A 3 / 2 fi W Z 3 . (33) 
3.5 First order global solution 

Let us remember the matching conditions we use in this paper: the metric 
components and their first derivatives have to be continuous through the 
hypersurface of zero pressure. Imposing these conditions on metrics fl2~T|) and 
fl3l!]) on the surface given by ([2~4l and (|25|) and keeping in mind the order of 
approximation we are concerned, an straightforward calculation leads to the 
following values for the exterior constants, 

A 2 w B 2 « 5 4 « , (34) 
these values for the interior constants, 

2^9 2 
m ^3, mjRi-1, j 1 &2+-n% j 3 » -- , 

a « 6 ~ «2 ~ &2 ~ (35) 

and this expression, cr ~ —5/6, for the unknown function in the surface 
equation. They all ensure that the matching has been properly done. 
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The above results lead to the following simple expression for the global 
metric: 



g int « (-1 + 3A - Xr, 2 ) T + (1 + 3A - Xrf) D - \Q 2 V 2 (T 2 + D 2 ) 



4 Second order metric in harmonic coordi- 
nates 

4.1 Non-linear exterior solution 

First of all we have to compute the right hand side of the equation ( fl~3l) . 
that is H and S. We just need approximate expressions of these quantities 
which can be obtained by means of the first order metric (1361) we worked out 
in the preceding section. In a standard post-Minkowskian approximation, it 
deserves nothing but a direct calculation keeping in mind we have only to 
retain terms of order less or equal to A 2 . However let us remember we have 
made a slow rotation assumption, so we also have to drop out terms that 
depend on powers of Q higher than Q 3 . Thus we get 





(36) 





x 2 n 2 



(V 3 -3W 3 ) , 



(37) 
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5r ?7 5 



where 



V n = P*(cos0) kidx { (n > 1) , 
W n = P n (cos6) edx 1 (n > 0) 



(38) 



are spherical harmonic vectors. 
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Substituting these expressions of S and H into equations ffTB"]) . we have 
a linear system of partial differential equations. It is a rather easy task to 
find a solution of this problem applying standard techniques, we choose that 

a 2 / 2 \ a 2 / n 2 \ 

^^- 2 ^( T °-3 Do )-3^( 1 - 2 vj F2 

X 2 Q 2 ( 32 1 1 \ 

+ — r- 2T 2 D 2 --£ 2 + — F 4 . 39 

r? 4 V 21 2 7 14 V 

This is not the best choice we can make because we need a lot of arbitrary 
constants if we attempt to match an exterior to an interior metric even in an 
approximate way. The linear exterior solution ( fl5l) . that is a solution of the 
homogeneous system ( TTBl) . give us all the constants we need. Therefore, our 
second order approximation to the exterior metric is a sum of three terms: the 
first order global exterior metric (IHUj) outside the source, a part of the solution 
to the exterior linear problem, mainly (l2"Tj) after having dropped out the 
Minkowski metric and the Z\ and Z% terms and having multiplied the result 
by A to obtain a second order expression, and the solution to the second order 
post-Minkowskian system (1391) . In fact, the first and second contributions to 
that expression can be obtained by substituting every constant in formula 
f[2"Tj) by its value after first order matching (|3~4"1) plus a new constant times 
A (we give to these new constant the same name since they can be seen as 
corrections to the previous ones); that is, 

Afo^l + AMo, M 2 ^-~ + AM 2 , 

A 2 ^XA 2 , B 2 ^XB 2 , B 4 ^XB 4 , (40) 

Ji and J 3 have not to be corrected because we stop our calculation at order 
A 2 . We arrive at the following expression for the exterior metric, 

~ ]rr o+ 1 + 2- + 2A 2 — )D 
V V V J 

D 2 ) + X 2 ^F 2 

X 2 Q 2 ^E 2 + X 2 Q 2 ^F i + h mh 
x 3 / 2 n /2 1 ,\ „ A 3 / 2 fi 
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+ 2'—^- [ - + -Q 2 ) Zi -2- — ^Z 3 . (41) 
T] 2 V5 3 / 7r] 4 s K 1 

As far as we have redefined the first order solution to introduce the second 
order solution constants, we can say that there are not other constants in our 



14 



approach than those of the linear solution. Moreover, they can be written as 
double series in powers of the source parameters A and Q. This is the kind of 
expressions we shall find for the multipole moments of the vacuum exterior 
metric. 



4.2 Non-linear interior solution 

First of all let us say something about the matching surface. The normaliza- 
tion factor (Tj0) up to order A 2 reads 



J fit + kit 



V^l + -A(/ tt + ftysin 2 #)+-A 

+ 2tlf tip V sm6 + Q 2 + / w ) V 2 sin 2 9 



(42) 



where k tt is the A 2 term in the expansion of 7^. Following what we did at 
first order, we use the continuity of the metric components on the matching 
surface to obtain an equation for this surface from the above expression. 
From the formula (HTj) we obtain the following expressions for the exterior 
metric components we need: 

2 n 2 

/«« gP 2 (cos0), 

k tt « - f M - - J + 2^r ( M 2 + - j P 2 (cos#) , 
V \ V/ TV V/ 

f tv * (cos 9), 

U^^-^P 2 (cos9). (43) 



We substitute them into ( 1421) and search for an explicit equation of the match- 
ing surface as a slight perturbation of the first order matching surface ff24|) . 



r 



1+ l-g+Ao-) £rP 2 (cosi 



(44) 



(the notation goes in the same way we have been using). Then, we find 

a w -\ + -M + M 2 . (45) 
5 o 

The value of ip on the matching surface up to the order we need is already 
known, we have just to set Mq = 1, its first order matching value, in its 
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expression f[2"""]) . Moreover if) is known too. We simply read out f tt form 
the expression of the first order global metric inside the source ( 1361) and 
substitute it into ( 12*31 to obtain it. Therefore we also have the expression 
for the pressure we need to write down the energy-momentum tensor up to 
order A 2 . That is, 

V « i ((3 - 2fi 2 )(l - V 2 ) - 5QyP 2 (cos0)) . (46) 



Let us proceed to obtain the second order interior solution. First we 
evaluate the right hand side of equation (1131) using the known first order 
solution ( l36l) . Now this task implies not only to obtain the non-linear terms 
in the metric but the energy-momentum tensor up to second order, i. e. 
neglecting terms in A 5//2 and f2 4 . The result is the following: 




o 



- (33 + 2Q 2 - fry 2 + ^W) D + 4^(1 - 4Q 2 ) V 2 F 2 
Tq \ 3 3/ 3ro 

\2q2 

+ — 5 -r7 2 (15T 2 + 23 J D 2 + 6^ 2 ) , 

H « 2—77 (3 - l& - v 2 + ^-Q 2 V 2 ) V x - 2- v (3 + 3Q 2 
r V 2 10 / r \ 

- r? 2 - l&rA W 1 - 2^ 3 (V 3 - 3 W 3 ) (47) 

Then, we build up the second order interior solution with the same ingredients 
we did the exterior solution: we obtain an exact solution of the inhomoge- 
neous system defined by equations (fl""|) and (T471) . 

^ inh = AV (-| + fi 2 + ^ 2 - IflV) T 

- a v (y + ^ 2 - h 2 + M D ° 

+ AV(4 + ^ + ^-^y)^ 2 

+ 3A W (-| + \ v 2 ) E 2 

+ ^A 2 OV ^15 T 2 + 23 D 2 + i F 4 ) , (48) 
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to which we add an adapted version of the general solution of the homoge- 
neous system and the first order interior solution ([3"B"]) . Thus we get 



g int « (-1 + 3A - Xr, 2 + A 2 m ) T + (1 + 3A - X V 2 + X 2 m ) D 
- Aft 2 (1 - Am 2 ) t] 2 (T 2 + D 2 ) + A 2 a £> + A 2 6 if 
+ A 2 ft 2 a 2 r/ 2 E* 2 + X 2 Q 2 b 2 r] 2 F* 2 + h inh 

+ 2X^ 2 n v ( i + -n 2 - - v 2 ) Zx - -A 3 / 2 fty z 3 . (49) 
v 3 5 / 7 

Here the free constants can also be thought as corrections of order A to the 
the first order interior solution constants, in the same way we suggested when 
we were solving the non-linear exterior problem, 

m — > 3 + Xm , m 2 — > — 1 + Am 2 , 

a — > Aa , &o — ► A6 , «2 — >■ Aa 2 , b 2 — > A6 2 , (50) 
ji and j3 do not need to be changed. 

4.3 Second order global solution 

Second order interior metric (1491) can be matched to the exterior metric (j4ip 
on the surface defined by (1441) and (1451) if an only if 

2 69 
M ^3 + -ft 2 , M 2 «-— , 
5 70 

^2~0, £ 2 « A (l + 1ft 2 ) , B A fH~ (51) 

and 



35 V 6 / 63 



m « i(21 + 2ft 2 ), m 2 ^-^j, 

a ^(21 + 2ft 2 ), 6 ~0, a 2 «-^, & 2 « . (52) 

3 15 105 

Substituting (|5*T1) into fj4"5"l) . we obtain the definite equation of the matching 
surface, 



l-^l-^A)o 2 A((os< 



The global second order solution reads 



(53) 



-1 + 3A - X V 2 + jX 2 (21 + 2ft 2 ) - A 2 (| - ft 2 ) 77 



2 



17 



!>2 

+ 2 A 



1 

2 
11 



- A — + -Q rj 



2 

2„2 



\ 2 2 

— A 77 



T + 

6 

1 H A Ar? 2 

70 14 7 



in 



7 



2 

21 



— 1 



1 + 3A - Ar/ 2 + 
- A^V 

An 2 



2 5 



7 

12 



A 2 (21 



2ft 2 



, 1079. 23 ^ 2 

1 H A Ar? 2 

210 14 ' 



r> 9 



1 AW (4 - 3^ 2 ) £ 2 + ^A 2 ft 

z 1 



n 



2 - 4 f 4 



7/ 



+ 2A 3 / 2 ^ f 1 + 1^ - Irf 
V 3 5 



;a 3 /W z 3 



9 



cxt 




69, A 

1 + — A + 32 

35 2I77 



£>9 



(54) 

Let us remark that the first order metric plays a relevant role in the above 
matching. When we evaluate the first order metric or their derivatives on the 
matching surface, substituting r by its expression (jSJ), there are generated 
terms of order higher than A that can not be ignored in the matching process. 

To end this subsection let us say something about the constants we have 
used to make the Lichnerowicz matching. The more interesting among them 
are obviously M , M 2 , J\ and J3. They are Thorne-Geroch-Hansen multipole 
moments, so they form part of an intrinsic (coordinate independent) char- 
acterization of vacuum asymptotically flat space-times. We must go back 
to the very beginning of the paper and try to join all the pieces into which 
we have split those constants if we want to write down them as they were 
introduced in formula (jT5|) . The tour starts in equation (12171) . goes through 
formulas (j3"4"|) . (j4"Ul) . and (loTj) to arrive at 

M « Ar 1 + A (3 + jjft 2 
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Jt « - 5 \ 3/2 nrl (l + , 

Js^-^n'rl (55) 

and 

A 2 ^0, B 2 « ^A 2 r 3 (l + V) , 5 4 ^-A A ^V 5 . (56) 

Even though A 2 , i? 2 and B4 are not intrinsic, they may be of some interest 
since they are set up in an unambiguous way by the Lichnerowicz matching 
conditions: they may carry some extra information about the source of the 
gravitational field, that is to say the kind of matter we are considering. 



4.4 Schwarzschild and Kerr metrics 

The static limit of our metric 



ft* - ("I + 3A - W + |A 2 - \\W + To 
+ (l + 3A - Xrj 2 + ^A 2 - H A y + IaV) D 

/ A A 2 A 2 \ / A A 2 A 2 \ 

g cxt ^ -1 + 2- + 6- -2- T + 1 + 2- + 6- + 4— ) D 

y r] 7] t] z 1 y 77 77 3?// 



A 2 /l 



r? 2 y 3 3577 



^2 , (57) 



must be an approximation to the global metric that results from matching 
Schwarzschild's metric to its well known interior metric (an exact solution of 
the Einstein equations for a perfect fluid having constant mass density) [24J . 
However, this metric is usually written in standard coordinates, therefore a 
direct comparison of this metric to (IBTj) is not allowed. Nevertheless, it has 
recently been shown [2B] that Schwarzschild's metric admits a global system 
of harmonic coordinates. Taking into account some formulas given in [25J, we 
can work out an approximate change of the radial coordinate which brings 
the Schwarzschild metric in standard coordinates to the form shown in ( IBTj) . 
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It reads, 



R(rj) 



r r] 



r r] 



10 



14 



A A 2 / 4 
! + " + - 3 + ^ 



V < 1 
T] > 1 



(5f 



where i? is the standard radial coordinate. 

Equation (1551) also relates the source radius in standard coordinates to 
the source radius in harmonic coordinates, that is, 



/ 1 0Q \ 

Ro = R(l)^r + \ + — A 2 J 



(59) 



We need this information to rewrite the Schwarzschild metric in harmonic co- 
ordinates, even outside the source, since the Schwarzschild mass is a function 
of Rq, 

(60) 



M = -iruRl « Ar (l + 3A) . 



M is equal to the mass monopole moment M as expected. 

We can establish that Kerr's metric does not fit our exterior metric by 
comparing the first multipole moments of the former (TSJ EE] , 



M = m , M 2 



-ma 



J\ = ma , J3 



-ma , 



(61) 



(here m and a stand for the usual parameters of the Kerr metric) to those of 
our metric. They cannot be matched for any choice of the free parameters 
A, Q and r : the following combination of multipole moments, 



A _ 3 A „ -I^/^ro 
M M 2 35 



(62) 



does not vanish in our approach but it is equal to zero if we evaluate it using 
the Kerr multipole moments. 



4.5 McLaurin ellipsoids 

The linear equation (|28|) that must be satisfied by the component of the 
metric deviation h u looks like the equation for the Newtonian gravitational 
field of a compact body with constant mass density \i. The classical solution 
to this problem is given by Poisson's integral; that is, a scalar field which 
tends to zero at infinity and such that the field and its first derivatives are 
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continuous on the surface of the body. Moreover, we know that the field 
outside of the body is uniquely determined by some constants, the multipole 
moments of the mass distribution. 

Let us consider an axisymmetric ellipsoid whose center is at the origin 
of the coordinate system and its axis lies on the z-axis. Then its mass and 
quadrupole moment are 

M = — /zerfc, 

Q = i M (b 2 - a 2 ) , (63) 
5 

where a and b stand, respectively, for the length of the ellipsoid semi-axis on 
the xy-plane and on the z-axis. These two parameters are constrained in a 
self-gravitating rotating ellipsoid by the equation (see for instance |27j ) 

5(1 + 25 2 ) 35 2 2 - 

(1-^)3/2 arccos 5 ~ ~ 3 n = ° ' (64) 

where 5 = b/a, and Q is related to the constant angular velocity of rotation 
of the ellipsoid u in the same way as the control parameter we have been 
using. Solving that equation for small values of Q, 

S^l-^n 2 , (65) 

and introducing for convenience a new parameter r$ defined by a = tqS^ 1 ^ 3 , 
we obtain 

= = °' ' 

Q = hr 3 5- 2 / 3 (5 2 - 1) w -hn 2 r 3 . (66) 

If we identify the parameters /x, us and r of the ellipsoid with those of our 
approximation scheme, these two expressions become identical with the A 
terms of the multipole moments Mq and Mi corresponding to our approxi- 
mate metric. Therefore, we can say that the Newtonian limit of the metric 
(|54p describes the gravitational field of a slowly rotating McLaurin's ellipsoid. 

Let us point out that the first dynamical multipole moment J\ is just the 
McLaurin ellipsoid angular momentum, 

J = 2 -Ma 2 uj = h 3 / 2 nr 2 5- 2 / 3 « ^A 3 / 2 fir 2 (l + ^ , (67) 

and J3 can also be related to the "quadrupole moment" of the angular mo- 
mentum density of the McLaurin ellipsoid. 
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5 Quo-harmonic coordinates 

Quo-harmonic coordinates were introduced some years ago by LI. Bel (see 
[26] for a brief review). They share with harmonic coordinates the property of 
being an admissible set of local coordinates for any space-time, since they are 
both defined by differential constraints that have sense whatever the metric 
is. However, quo-harmonic coordinates unlike harmonic coordinates require a 
timelike vector field to be properly defined. That field is used to set the time 
coordinate and to introduce the projector into the plane orthogonal to the 
field (quotient metric), which in turn provides spatial coordinates as three 
independent harmonic functions of that three dimensional positive definite 
Riemannian metric (quo-harmonic coordinates). If these functions are chosen 
appropritely, they may look like Cartesian coordinates in the same grounds 
as harmonic coordinates can be thought of as a generalization of Minkowski's 
coordinates to curved space-times. 

It is sometimes possible to choose these coordinates as first integrals of 
the timelike vector field. That property singles out a class of vector fields 
in any Riemannian manifold [28]. Those fields yield to introduce what Bel 
calls quo-harmonic congruences, an attempt to understand rigid motion in a 
framework less restrictive than that of Born's [291 SO] • 

5.1 Quo- harmonic coordinates and post-Minkowskian 
approximation 

The timelike Killing vector £ defines a quo-harmonic congruence in any space- 
time whose metric is given by (CD). The assumptions we have made on the 
coordinates (see section 2) identify the t coordinate with the time we can 
associate with that Killing field, since £ has canonical form in this system of 
coordinates. This coordinate is an harmonic function of the four dimensional 
metric. This constraint on the time coordinate is not required by the quo- 
harmonic scenario, which does not define this coordinate, but it is given as a 
bonus by the Papapetrou structure we assume for the metric. Moreover, the 
projector operator associated with £ 



does not depend on t: it is a true three-dimensional metric. This is more 
evident in our notation, 



9t = 9-g(€,€) 





(69) 
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where {x 1 } are the Cartesian-like coordinates we introduce in section 2. Then, 
if we want them to be quo-harmonic coordinates, we have to impose the 
following condition on the Riemannian connection of the three dimensional 
metric ( 1691) 

9^k = (70) 
From the point of view of a post-Minkowskian perturbation approach, 
this coordinate condition is as admissible as the harmonic condition ( JTU1) . We 
obtain a well defined system of differential equations that can be formally 
solved by iteration. The first-order approximation to the metric unambigu- 
ously determines the inhomogeneous part of the system of equations which 
can be solved to work out the second-order metric and so on. 

The starting point is a system of linear partial differential equations that 
are slightly more complex than those mentioned before, say (TH|) or (1281) . 

A/ioo = -167rfoo , 

Ah 0i = -167r£ j , 

Ahij - dijhoo = -16-Ktij , 

d k (h kt -^5 kt 6 lm h lm ) =0. (71) 

Their solution can be written using the basis of spherical harmonic tensors 
we introduced in the harmonic case. Outside the source, t a p = 0, it reads, 

~ M n ( n + 2 1 \ 

h = 2 Y — ^ T n + 2 D n + — r F n+2 

^ r"+! \ n 2n + 3 2(2n + 3) + / 

OO T OO A OO 75 

n=l ' n=l ' n=2 ' 

inside the source, it reads 

(2 1 
h = h inh + m T + m l rT 1 + m 2 r 2 f T 2 + - D 2 - - H 

+ m 3 r 3 (t z + ^D z -^{h x -H\ 

°° / 77 — 1 1 \ oo 

+ Y m n r n T n + 2 D n - F* 2 + Y Jn r n Z n 

^ 4 \ 2n-l 2(2n-l) n - 2 J^^ Jn 



+ a D + b H + Y a n r n E* n + Y b n r n F* n . (73) 

n=l n=2 

The notation for the constants is the same as that used in the harmonic 
case because they are required to have the same structure on A and Q we 
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assumed there: constants that are not compatible with equatorial symmetry 
must be set equal to zero, the other constants are functions of the physical 
parameters, A, fl, and ro like those we introduced by formulas ( !20l) and ( j32l ). 
This implies that t is given by ( 1271) up to the first order in A (even to order 
A 3 / 2 ). Though its formal expression coincides with that of the harmonic case, 
of course it does not lead to the same inhomogeneous solution, 



h 



inh 



Dr 



15 



6 



(74) 



Following the scheme we introduced in the harmonic case, we obtain the 
first-order exterior and interior metrics by substituting ( 1201) and ( 1321) into 
expressions ( 1721) and ( 1731) and neglecting terms in A 2 and Q 4 , we find out the 
equation of the matching surface at this order of approximation, we impose 
the continuity of the metric on that surface, etc. . . The procedure is as before. 
Therefore, we shall merely report the final results and comment on them. 



5.2 



Second-order global solution 

The metric can be written in quo-harmonic coordinates as follows: 
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(75) 



The equation of the matching surface, on which the metric and their first 
derivatives are continuous up to the order of approximation we are concerned, 
is 



r 



22 
35' 



1 - — A ) ft 2 P 2 (cos( 



(76) 



The pressure of the fluid has the same expression we found in the harmonic 
case but obviously the meaning of the coordinates and the parameters is not 
the same. It suggests that quo-harmonic coordinates differ from harmonic 
coordinates at least in first order terms. 



5.3 Multipole moments 

Quo-harmonic coordinates are asymptotically Cartesian and mass-centered 
coordinates in the sense of Thorne [TS] (see Appendix). Then we can read 
the multipole moments of the vacuum exterior field from the g a and g tip 
components of the metric ( J751 . However we do not need to do this search 
because they must be equal to the constants of the exterior metric that 
appear in these two components. Lichnerowicz's matching provides for them 
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the following values: 



M « Ar c 



M 2 « --AflVg ( 1 



1 + - (21 + 2Vt 
5 V 

139, 



35 



J X » ^A 3 / 2 fir 2 (l + j^ 2 ) 



J 3 



A 3 /W 



(77) 



These expressions look like those we obtained for the multipole moments 
in the harmonic case but differ from them in A 2 terms. Nevertheless, mul- 
tipole moments are coordinate independent quantities so both expressions 
must lead to the same numerical values for any set of the source parameters 
we choose. Let us be more precise. The parameters /i and u are clearly 
constants that characterize the fluid but ro, the mean source radius, is a 
constant of a different class: it can be defined as the value taken by the 
radial coordinate on the matching surface corresponding to the static limit 
of our metric. Therefore r does not have to be the same in quo-harmonic 
coordinates and in harmonic coordinates. However they both must be linked 
if they parametrize the same surface of the space-time; that is, if we are 
talking about just a single member of this one-parameter family of metrics. 
We can obtain some insight into the relationship between quo-harmonic and 
harmonic values of r which lead to the same metric by identifying the ex- 
pressions of multipole moments in both systems of coordinates. Equation 
(I77p can be transformed into (j55p by means of the following substitutions 



r 



roil 



-A 



A 



xn-ix 



(78) 



Let us note that A is a function of ro, then it has different values in one or 
another system of coordinates this is the meaning of the second substitution. 
On the other hand Q is just a function of fi and u so it does not change. 

The above comment points to the change from quo-harmonic to harmonic 
coordinates for the Schwarzschild metric to be in the origin of the substitution 
rule flTHj) . That is true. If we substitute the quo-harmonic radial coordinate 
in the static limit of the quo-harmonic metric (|75p up to order A, 



g int « (-1 + 3A - Arr 2 ) T + ( 1 + 4A - ^X V 2 ) D ~A/f F 2 . 
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■1 + 2- To + 1 



5 377J 0+ 3?7 l 1 




(79) 
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in the following way, 



r 1 

2 




) 



for C/int 



r 



(80) 



r 1 — A — 

2r 




) 



for gf, 



ext 



which implies that r$ and A have to be substituted at the same time in the 
way we suggested in (ITS]) . If so we reach the static limit for the harmonic 
case metric fl57|) . 

The substitutions (|78p make no sense when we apply them to other exte- 
rior or interior constants. In fact, A 2 , B 2 , B 4 , a , bo, a 2 and b 2 , which were 
at least of order A 2 in the harmonic of order A in the quo-harmonic 

case. It makes all the calculations more cumbersome and the final expression 
of the metric more involved, as one can see by comparing ( 1751) to ( 15^1) . 
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Appendix: Approximate metric 

Here we write down the components of the metric in the original Papapetrou 
form (00). This is a standard way to introduce stationary axisymmetric met- 
rics. Therefore, some interesting properties of the metric structure can be 
easier recognized looking at the following expressions than they could be by 
examining the formulas of the preceding sections, even though they are more 
suitable to run the perturbation algorithm. We use the obvious shorthand 
notation P t m = Pj m (cos0). 

Harmonic coordinates 

A.l Interior metric components in harmonic coordinates 



ltt = -l + 3A + -A 2 (21 + 2Q 2 
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A. 2 Exterior metric components in harmonic coordinates 
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Quo-harmonic coordinates 

A. 3 Interior metric components in quo-harmonic coordinates 
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